Maximally entangled mixed states of two qubits 
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We consider mixed states of two qubits and show under which global unitary operations their 
entanglement is maximized. This leads to a class of states that is a generalization of the Bell states. 
Three measures of entanglement are considered: entanglement of formation, negativity and relative 
entropy of entanglement. Surprisingly all states that maximize one measure also maximize the 
others. We will give a complete characterization of these generalized Bell states and prove that 
these states for fixed eigenvalues are all equivalent under local unitary transformations. We will 
furthermore characterize aU nearly entangled states closest to the maximally mixed state and derive 
a new lower bound on the volume of separable mixed states. 
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In this paper we investigate how much entanglement 
in a mixed two qubit system can be created by global 
unitary transformations. The class of states for which 
no more entanglement can be created by global unitary 
operations is clearly a generalization of the class of Bell 
states to mixed states, and gives strict bounds on how 
the mixedness of a state limits its entanglement. This 
question is of considerable interest as entanglement is 
the magic ingredient of quantum information theory and 
experiments always deal with mixed states. Recently, 
Ishizaka and Hiroshima [Q independently considered the 
same question. They proposed a class of states and con- 
jectured that the entanglement of formation |^ and the 
negativity ||] of these states could not be increased by 
any global unitary operation. Here we rigorously prove 
their conjecture and furthermore prove that the states 
they proposed are the only ones having the property of 
maximal entanglement. 

Closely related to the issue of generalized Bell states 
is the question of characterizing the set of separable den- 
sity matrices [||, as the entangled states closest to the 
maximally mixed state necessarily have to belong to the 
proposed class of maximal entangled mixed states. We 
can thus give a complete characterization of all nearly 
entangled states lying on the boundary of the sphere of 
separable states surrounding the maximally mixed state. 
As a byproduct this gives an alternative derivation of the 
well known result of Zyczkowski et al. Q that all states 
for which the inequality Tr(yO^) < 1/3 holds are separa- 
ble. 

The original motivation of this paper was the follow- 
ing question: given a single quantum mechanical system 
consisting of two unentangled spin-1/2 systems, i.e. two 
qubits, one in a pure state and another in a maximally 
mixed state, does there exist a global unitary transfor- 
mation on both qubits such that they become entangled? 
Surprisingly, the answer is yes. In this paper we solve the 
more general question: how can one maximize the entan- 
glement of an arbitrary mixerd state of two qubits using 
only unitary operations. If not only unitary operations 
but also measurements were allowed, it is clear that a Von 



Neumann measurement in the Bell basis would immedi- 
atly yield a singlet. Here however we restrict ourselves to 
unitary operations. Obviously, these unitary operations 
must be global ones, that is, acting on the system as 
a whole, since any reasonable measure of entanglement 
must be invariant under local unitary operations, acting 
only on single qubits. As measures of entanglement, the 
entanglement of formation (EoF) the negativity j|] 
and the relative entropy of entanglement [Q were chosen. 

The entanglement of formation of mixed states is de- 
fined variationally as Ef{p) = minj^.j 'YliiPi^ii^i) where 
p — ^iPiipiipi- For 2x2 systems the EoF is well- 
characterized by introducing the concurrence C : 



Ef{p)^f{C{p)) = H 



1 + \/l - C2 



C{p) = max(0, CTi - (72 - - 0-4). 



(1) 
(2) 



Here {17^} are the square roots of the eigenvalues of the 
matrix A arranged in decreasing order 



A — pSp* S 



(3) 
(4) 



H{x) is Shannon's entropy function and ay is the Pauli 
matrix. It can be shown that /(C) is convex and 
monotonously increasing. Using some elementary linear 
algebra it is furthermore easy to prove that the numbers 
{(7i\ are equal to the singular values |^ of the matrix 
S^/p. Here we use the notation y/p = $A^/^ given 
^A^^, the eigenvalue decomposition of p. 

The concept of negativity of a state is closely related 
to the well-known Peres condition for separability of a 
state . If a state is separable (disentangled) , then the 
partial transpose of the state is again a valid state, i.e. 
it is positive. For 2x2 systems, this condition is also 
sufhcient 0|. It turns out that the partial transpose of 
a non-separable state has one negative eigenvalue. From 
this, a measure for entanglement follows: the negativity 
of a state H is equal to the trace norm of its partial 
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transpose. We will adopt the definition of negativity as 
twice the absolute value of this negative eigenvalue: 



En{p) = 2max(0, — A4), 



(5) 



where A4 is the minimal eigenvalue of p^^ . In the case 
of two qubits, this is equivalent to the trace norm of the 
partial transpose up to an affine mapping. 

The relative entropy of entanglement was proposed by 
Vedral and Plenio [Q as a measure of entanglement mo- 
tivated by the classical concept of KuUback-Leibler dis- 
tance between probability distributions. This measure 
has very nice properties such as being a good upper 
bound for the entanglement of distillation. It is varia- 
tionally defined as 



Er{p) 



min Tr {p log p — p log a) 



(6) 



where D represents the convex set of all separable density 
operators. 

We now state our main result: 

Theorem 1 Let the eigenvalue decomposition of p be 
p = <I>A<&^ 

where the eigenvalues {Xi} are sorted in non- ascending 
order. The entanglement of formation is maximized if 
and only if a global unitary transformation of the form 



U ^{Ui® U2) 



/ 1 \ 

I/V2 l/\/2 

1/V2 -1/V2 
Vox 0/ 



is applied to the system, where Ui and U2 are local uni- 
tary operations and is a unitary diagonal matrix. 
This same global unitary transformation is the unique 
transformation maximizing the negativity and the relative 
entropy of entanglement. The entanglement of formation 
and negativity of the new state p' = UpW are then given 

by 



Ef{p') = f (max(o, Ai - A 
Em{p) = max 



2V/A2X 



(0, v/(Ai - A3)2 + (A2 - XiY -X2- A4) 



respectively, while the expression for the relative entropy 
of entanglement is given by 

En{p') = Tr(plogp) - Ai log((l - a)/2) - 
A2 log((a + 6 + 2(A2-A4))/4)- 
A3 log((l - 6)/2) - A4 log((a + b - 2(A2 - A4))/4) 



a^{d- ^d^ - 4(1 - Ai)(l - A3)(A2 - A4)2)/(2(l - A3)) 
b={d- ^d^ - 4(1 - Ai)(l - A3)(A2 - A4)2)/(2(l - Ai)) 
d = A2 + A4 + (A2 - A4)^ 



The class of generalized Bell states is defined as the states 
p' thus obtained. These states are the maximally entan- 
gled mixed states (MEMS-states). 

We now present the complete proof of this Theorem. 
The cases of entanglement of formation, negativity and 
relative entropy of entanglement will be treated indepen- 
dently. We start with the entanglement of formation. 

As the function f{x) is monotonously increasing, max- 
imizing the EoF is equivalent to maximizing the concur- 
rence. The problem is now reduced to finding: 

Cmax = max (0, (Ti - (72 - 0-3 - 0-4) (7) 

(76(7(4) 

with {ai} the singular values of 

Q ^ Ai/2$T[/T5'C/$Ai/^ (8) 

Now, $, U and S are unitary, and so is any product of 
them. It then follows that 



Cmax < max (0, 0-1 - (72 - 0-3 - (74) 

V<£U(i) 



(9) 



with {ffi} the singular values of A^^'^VA^/^ . The inequal- 
ity becomes an equality if there is a unitary matrix U 
such that the optimal V can be written as ^^U'^ SU^. 
A necessary and sufficient condition for this is that the 
optimal V be symmetric (V — V'^): as S is symmetric 
and unitary, it can be written as a product Sf Si, with Si 
again unitary. This is known as the Takagi factorization 
of 5* HJ. This factorization is not unique: left-multiplying 
Si with a complex orthogonal matrix O {O^O — 1) also 
yields a valid Takagi factor. An explicit form of 5*1 is 
given by: 



^1 




(10) 



If V is symmetric it can also be factorized like this: 
V = ViVi. It is now easy to see that any U of the 
form 



U = S\OVi^\ 



(11) 



with O real orthogonal, indeed yields V = ViVi. 

To proceed, we need two inequalities concerning singu- 
lar values of matrix products. Henceforth, singular val- 
ues, as well as eigenvalues will be sorted in non- ascending 
order. The following inequality for singular values is well- 
known 01: 



Lemma 1 Let A M„,r(C), B e Mr,„(C). Then, 



Y^aUB)<Y^a^{A)am. 

i=\ i=l 

for k = 1, . . . ,q = min{n, r, m}. 



(12) 
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Less known is the following result by Wang and Xi : 

Lemma 2 Let A £ M„{C), B € M„^„(C), and 1 < ii < 
• • • < jfc < n. Then 



(13) 



t=i 



Set n = 4 in both inequalities. Then put A: = 1 in the 
first, and fc = 3, ii = 2, Z2 = 3,i3 = 4 in the second. 
Subtracting the inequalities then gives: 

aMB) - {(T2{AB) + (73 (AS) + (Ji{AB)) < 
ai{A)ai(B) - a2{A)a4{B) - g^{A)g^{B) - G4{A)g2{B) 

Furthermore, let A = A^/^ and B = Vh^l"^, with A pos- 
itive diagonal and with the diagonal elements sorted in 
non-ascending order. Thus, Gi{A) = Gi{B) = y/Xi. This 
gives: 

(ai-(a2 + (T3 + a4))(Ai/VAi/2) < Ai - (2 A3). 

It is easy to see that this inequality becomes an equality 
iff V is equal to the permutation matrix 



/ 1 

1 

10 

VO 1 0. 



(14) 



multiplied by an arbitrary unitary diagonal matrix D^. 
Therefore, we have proven: 

maxvec/(4) {(^1 - (^2 + <J3 + G4)){A^^^VA^/^) = 

Ai- (2VA^+A3). (15) 

We can directly apply this to the problem at hand. The 
optimal V is indeed symmetric, so that it can be decom- 
posed as V = V^V\. A possible Takagi factor is: 



/ 1 \ 
I/V2 

10 

\ -VV2 / 



(16) 



The optimal unitary operations V are thus all of the form: 
\J = SlOViD^J^^^ with O an arbitrary orthogonal ma- 
trix. It has to be emphasized that the diagonal matrix 
will not have any effect on the state p' = J7$A$t[/t^ 
To proceed we exploit a well-known accident in Lie 
group theory : 



SU{2)®SU{2) = 50(4). 



(17) 



It now happens that the unitary matrix Si is exactly 
of the form for making Si{Ui U2)Sl real for arbi- 
trary {Ui,U2} e SU{2). It follows that Si{Ui ® U2)Sl 



is orthogonal and thus is an element of 50(4). Con- 
versely, each element Q E SO (A) can be written as 
Q — Si{Ui (8) U2)Sl. On the other hand the orthogonal 
matrices with determinant equal to —1 can all be written 
as an orthogonal matrix with determinant 1 multiplied 
by a fixed matrix of determinant —1. Some calculations 
reveal that 



Si 



/ 1 
1 
1 

Vo 







-1 



)S\Vi 



/ 1 Q\ 
10 
0-10 

Vo 1/ 



We conclude that for each O G 0(4) and unitary di- 
agonal, there exist t/i, C/2 € SU{2) and I?^' unitary diag- 
onal, such that U = ^jOViL'^^t ^ {Ui®U2)SlViD^,<i>''. 

It is now easy to check that a unitary transformation 
produces maximal entanglement of formation if and only 
if it is of the form 



(Ui (g, U2) 



/ 

1/V2 
1/V2 
V 1 



1/V2 
-1/V2 

/ 



(18) 



This completes the proof of the first part of the Theorem. 

We now proceed to prove the second part of the Theo- 
rem concerning the negativity. This proof is based on the 
Rayleigh-Ritz variational characterization of the minimal 
eigenvalue of a Hermitian matrix: 



Amm (p 



'^'*^= min p'^^\x){x\ 

2;: I |a; 1 1 —1 

— mill p(\x){x\)^^ 

a^: I |a; 1 1 — 1 



(19) 



The eigenvalue decomposition of {\x){x\f'^ can best be 
deduced from its singular value decomposition. Let x de- 
note a reshaping of the vector a; to a 2 x 2 matrix with 



as follows: 



y • 



e^\x). Introducing the permutation matrix 
the partial transpose can be written 



{\x){x\) 



Ta 



(20) 



The proof of this statement is elementary. We denote the 
Schmidt decomposition of the vector \x) by 



(21) 



where the diagonal elements of S are given by g\^gi. 
Since x is normalized we can parameterize these as 
cos(a), sin(a) with < a < 7r/4 (to maintain the or- 
dering). We get 



{\x){x\) 



Ta 



Po{Ui(g)U2){^(E)^){U2(E)Ui 



(22) 



This clearly is a singular value decomposition. The ex- 
plicit eigenvalue decomposition can now be calculated us- 
ing the basic property of Pq that Po{A(g) B) = {B(g)A)Po 
for arbitrary A, B. It is then easy to check that the eigen- 
value decomposition of {\x){x\)'^^ is given by: 
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{\x){x\f^ =V{x)D{a{x))V{xy 



(23) 



where D{a{x)) is the diagonal matrix with eigenvalues 
(erf , (Ticr2, erf, -(Ti(T2) and 



V{x) = {Ui{x)®U2{x)) 



/ 1 \ 

1/72 l/\/2 

1/72 -1/V2 

yo 1 / 



(24) 



For the problem at hand, we have to minimize the min- 
imal eigenvalue of (UpW)'^^ over all possible U E C^(4). 
Thus, we have to minimize: 

minTr U'i>A<S>^U''V{x)D{a{x))V{x)^ 

U.x 

— min min TTAW^D{a)W, (25) 

a W 

where we have absorbed the eigenvector matrix $ of p, 
as well as V{x)'', into U, yielding W. Now, the min- 
imization over W can be done by writing the trace in 
components 



7(a) = Tr AW^D{a)W 
= J2d,{a)\W,,\'X. 



d{afj{W)X, 



(26) 



where d{a) and A denote the vectors containing the di- 
agonal elements of D{a) and A, respectively. J(W) is a 
doubly stochastic matrix formed from W by taking the 
modulus squared of every element. The minimum over 
all W is attained when J{W) is a permutation matrix; 
this follows from Birkhoff 's theorem j|] , which says that 
the set of doubly-stochastic matrices is the convex clo- 
sure of the set of permutation matrices, and also of the 
fact that our object function is linear. Since the compo- 
nents of (7 and A are sorted in descending order and A is 
positive, the permutation matrix yielding the minimum 
for any a is the matrix 



Jo 



/ 1 

10 

10 

V 1 . 



(27) 



Thus W has to be chosen equal to Jq multiplied by a di- 
agonal unitary matrix D^. Hence, the minimum over W 
is given by J2j=i Aj(i4+i_j(a). Minimizing over a gives, 
after a few basic calculations: 



cos(2a) 



Aa — A4 



v/(Ai - A3)2 + (A2 - A4)2 
9{a) = (a2 + A4 - v/(Ai-A3)2 + (A2-A4)2) /2. 

This immediately yields the conjectured formula for the 
optimal negativity. 



We now have to find the U for which this optimum is 
reached. As V^xYU^ = W, it follows that the optimal 
unitary transformation U is given by U = V{x)JoD^^'': 



( 0, 







C/ = (C/i ® U2) 



1\ 



l/\/2 l/%/2 
l/\/2 -I/V2 

V 1 / 



D^^^ (28) 



This is exactly the same U as in the case of entanglement 
of formation. 

Next we move to the third part of the theorem con- 
cerning the relative entropy of entanglement. We first 
prove two lemmas. 

Lemma 3 Consider the class of superoperators 
r(p)=^a,[/,pt//, 

i 

where all Ui are unitary, and the ai form a distribution. 
Then, for any state p that is invariant under T , we have 
for the relative entropy: 

S{p\\a)>S{p\\T\a)). 

Proof. The proof of this lemma is heavily inspired by 
theorem 6 in Rains Q. From S{p) = S{T{p)), we find 

5(p||a) =Trplogp-Trploga 

= Trplogp-TrT(/9) log a 

= Tr p log p - ^ a, Tr U^pUj log a 

i 

= Tr p log p - ^ a, Tr p log{u}aU - i) 

i 

> Tr p\og p ~ Ti plogiY, a^U| aU -i) 

i 

= Sip\\TH<y)), 

where in the penultimate line, we have used the subaddi- 
tivity of the relative entropy w.r.t. its second argument. 

Lemma 4 For p of the form p = UAU'^ with 



U = 



/ 1 \ 

1/72 1/72 
1/72 -1/72 

Vol / 



and A containing the ordered eigenvalues of p, 

Er{p) ^ min S'(pIIo-). 

aennMEMS 

where MEMS is the class of maximally entangled mixed 
states. 
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Proof. Define the following superoperator: 

T{p) = V di&g{V^ pV)V^ . 

Here, diag(p) is the superoperator that sets all ofF- 
diagonal elements of p equal to zero while keeping the 
diagonal ones intact. This superoperator can also be 
written as 

diag(p)==^P,pP,/2", 

i 

where Pi runs through all possible diagonal matrices hav- 
ing only +1 or —1 on their diagonal It follows that T 
is of the form mentioned in the first Lemma and, further- 
more, that it is a self-dual superoperator, i.e. = T. 

It is obvious that aU MEMS states (with Ui = U2 = 1) 
are left invariant by T. We will now show that any such 
T maps separable states to separable states. Consider 
thereto the pure product states only; if the proposition 
is valid for pure product states, it will be valid for all 
separable states (by linearity). The most general pure 
product state has the state vector = [ac,ad,bc,bd), 
with a, 6, c, d complex numbers. Then, since 

yt^ = ((arf 4_ &c)/%/2, bd, {ad - 6c)/V2, ac), 

r(V'V'^) = V Aisig{vHilj^V)V^ 
= VKV\ 

where the diagonal elements of A are, in order, 

{\{ad + bc)/V2\^, |6d|2, |(ad- 6c)/V2|2, \ac\^). 

As these values are not necessarily sorted, T{ipip'') need 
not be MEMS. However, it is still possible to apply the 
formula for the negativity of MEMS states which says 
that 



En{p) - max(0, V i^i - + i>^2 - X^)^ - A2 - A4). 

As can be easily checked, the validity of this formula does 
not rely on the ordering of the Ai, as long as each Ai per- 
tains to the i-th column of V. In particular, using 

Ai - A3 = 2n{ab{cd)*) 
A2±A4 = |M|2±|acp, 

we get for the negativity of T{il;ip''): 

EN{r{i)^'^)) = max(0,F) 

with 



F = v/(Ai - A3)2 + (A2 - A4)2 - (A2 + A4) 
= ^Jm{ab{cd)*Y + |6d|4 + \ac\^ - 2\abcdY 

= y/\bd\^ + |ac|4 + 2|a6cd|2 - 43(a6(cd)*)2 

= v/(|MP + |ac|2)2 - 43(a5(cd)*)2 - (|6d|2 + jacp) 
< 0. 



Hence, T{tj)ijj'^) is separable, as we set out to prove, so 
that T maps separable states to separable states. 

From the previous discussion it also follows that states 
of the form T^Ay^ are separable if and only if the eigen- 
values satisfy 



V(Al - A3)2 + (A2 - A4)2 - A2 - A4 < 



(29) 



Furthermore, states l^AT^^ are obviously invariant under 
T . Hence, letting a traverse all separable states of this 
form generates the same set of states T(cr) as letting a 
traverse all separable states without restriction. There- 
fore, 



Eb,{p) = minS'(pllcr) 



>min5(H|T(a)) 



min S'(p||cr) 

--VKVeD 



Comparing the first and the third line, we immediately 
see that the inequality must be an equality. 

Actually, an even stronger result holds, as we can 
restrict ourselves in this minimization to states a = 
Vhy"^ G D where the diagonal elements appear in de- 
scending order (Ai > A2 > A3 > A4). In other words, 
a may be taken from the set of separable MEMS states. 
To see this, note that, as p and a are both MEMS, 

S{p\W) = ^p,;(lgpj - IgAi), 



where the pi are the sorted eigenvalues of p, and Ai are 
the not necessarily sorted eigenvalues of cr. It is easy to 
see that one always gets a lower relative entropy by per- 
muting the Ai into descending order This ends the 
proof of the lemma. 

It is now easy to prove the last part of the main the- 
orem: Because the a are restricted to separable MEMS 
states, this means that, for any global unitary U, Uall'' 
is still separable. Hence, for p G MEMS, 

Er(p) = min S{p\\a) 
^ ' .GDnMEMS 

min SiUpU^WUdU^) 
^ecnMEMS 

> min SiUpU^WUaU^) 
= Er{UpU^), 

where the inequality in the penultimate line arises 
because the minimization domain has been enlarged. 
Therefore the MEMS states have larger relative entropy 
of entanglement then all states that can be obtained from 
it by doing global unitary operations. 
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The explicit calculation of the relative entropy of en- 
tanglement of the maximally entangled mixed states is 
now a tedious but straightforward exercise, whose result 
is quoted in the theorem. This completes the proof of 
Theorem 1. □ 

Let us now analyze more closely the newly defined class 
of generalized Bell states. We already know that U is 
unique up to local unitary transformations. It is easy 
to check that the ordered eigenvalues of the generalized 
Bell states for given entanglement of formation /(C) are 
parameterized by two independent variables a and /3: 

< a < 1 




a, ^/S- 



V2a) 



Ai = 1 



Ao = 



A., = 



l-C 



6 




1 - 


C 


6 




1 - 


C 


6 



[a + %/2/3)2 
(3- (V2a + /3)2) 



(30) 



For given EoF there is thus, up to local unitary transfor- 
mations, a two dimensional manifold of maximally entan- 
gled states. In the case of concurrence C = 1 the upper 
and lower bounds on (3 become equal and the unique 
pure Bell states arise. Another observation is the fact 
that A4 of all generalized Bell states is smaller then 1/6. 
This implies that if the smallest eigenvalue of whatever 
two-qubit state exceeds 1/6, this state is separable. 

A natural question is now how to characterize the en- 
tangled states closest to the maximally mixed state. A 
sensible metric is given by the Frobenius norm 1||2 = 
'^i ~ 1/4- This norm is only dependent on the 
eigenvalues of p and it is thus sufficient to consider 
the generalized Bell states at the boundary of entangled 
states where both the concurrence and the negativity be- 
come zero. This can be solved using the method of La- 
grange multipliers. A straightforward calculation leads 
to a one-parameter family of solutions: 

< x < i 




A, = :r 



(31) 



The Frobenius norm — 1||2 for all these states on the 
boundary of the sphere of separable states is given by 
the number \J\jVl. This criterion is exactly equiva- 
lent to the well-known criterion of Zyczkowski et al. : 



Trp^ = 1/3. Here, however, we have the additional 
benefit of knowing exactly all the entangled states on 
this boundary as these are the generalized Bell states 
with eigenvalues given by the previous formula. Further- 
more, Zyczkowski et al. |^ proposed a lower bound on 
the volume of separable states by considering the ball 
of states that remain separable under all global unitary 
transformations. Clearly the criterion Af < 1/3 can 
be strengthened to Ai — A3 — 2\J\2\a < 0. Some tedious 
integration then leads to a better lower bound for the 
volume of separable states relative to the volume of all 
states: 0.3270 (as opposed to 0.3023 of @). 

Further interesting properties of the maximally entan- 
gled mixed states include the fact that the states with 
maximal entropy for given entanglement all belong to 
this class. This can be seen as follows: the global entropy 
of a state is a function of the eigenvalues of the density 
matrix only. Therefore the states with maximal entangle- 
ment for given entropy can be found by first looking for 
the states with maximal entanglement for fixed eigenval- 
ues, followed by maximizing the entropy of the obtained 
class of (maximally entangled) mixed states. 

In conclusion, we have generalized the concept of pure 
Bell states to mixed states of two qubits. We have proven 
that the entanglement of formation, the negativity and 
the relative entropy of entanglement of these generalized 
Bell states could not be increased by applying any global 
unitary transformation. Whether their entanglement of 
distillation is also maximal is an interesting open prob- 
lem. 

We thank Lieven Vandersypen for bringing the prob- 
lem under our attention, and J. Dehaene, L. De Lath- 
auwer and K. Zyczkowski for valuable comments. T. De 
Bie is a Research Assistant with the Fund for Scientific 
Research-Flanders (FWO- Vlaanderen) . 
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